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In this article we study two problems that arise when using timing and amplitude estimates 
from a network of interferometers (IFOs) to evaluate the direction of an incident gravitational wave 
burst (GWB). First, we discuss an angular bias in the least squares timing-based approach that 
becomes increasingly relevant for moderate to low signal-to-noise ratios. We show how estimates of 
the arrival time uncertainties in each detector can be used to correct this bias. We also introduce 
a stand alone parameter estimation algorithm that can improve the arrival time estimation and 
provide root-sum-squared strain amplitude (hrss) values for each site. In the second part of the 
paper we discuss how to resolve the directional ambiguity that arises from observations in three 
non co-located interferometers between the true source location and its mirror image across the 
plane containing the detectors. We introduce a new, exact relationship among the hrss values at 
the three sites that, for sufficiently large signal amplitudes, determines the true source direction 
regardless of whether or not the signal is linearly polarized. Both the algorithm estimating arrival 
times, arrival time uncertainties, and hrss values and the directional follow-up can be applied to any 
set of gravitational wave candidates observed in a network of three non co-located interferometers. 
As a case study we test the methods on simulated waveforms embedded in simulations of the noise 
of the LIGO and Virgo detectors at design sensitivity. 

PACS numbers: 



I. Introduction 

As the first generation of gravitational wave interfer- 
ometers perform observations at or near their design sen- 
sitivities, new methods are being developed to detect and 
characterize gravitational wave burst signals. Accurate 
determination of the source direction is fundamental for 
all analyses of a candidate signal. 

There are two approaches to localizing a source on the 
celestial sphere: coincident and coherent. Coincident 
methods [l|] analyze the data from each detector sepa- 
rately and then identify events that occur simultaneously 
in multiple interferometers. The arrival times and am- 
plitudes in each detector can then be used to determine 
the source direction if the signal is linearly polarized or 
has polarization peaks separated by less than the timing 
uncertainty. Coherent methods 0, d, |3j IE] combine the 
data streams from multiple detectors into a single statis- 
tic before reconstructing gravitational wave candidates. 

One previous coincident study on source location esti- 
mation was completed by Cavalier et al[l[. In that work, 
arrival times (assumed known) and their uncertainties 
(assumed gaussian) were input to a minimimization 
routine to determine the source location. The approach is 
applicable to an arbitrary number of interferometers and 
shows excellent resolution in Monte Carlo simulations 
with arrival time uncertainties on the order of 0.1 ms. 
In this paper we implement an equivalent least squares 
approach in the two and three interferometer cases, but 
also consider arrival time uncertainites up to 3 ms. These 
larger uncertainties have been observed in studies with 
real noise@, 0| and reveal a systematic bias in the source 



directions obtained from timing reconstruction methods 
for moderate to low signal-to- noise ratios (SNRs). We 
define the SNR p as 

where 



h{f) = / h{t)e-'^''f'dt (2) 

is the Fourier Transform of the gravitational wave sig- 
nal and 

S{f)^Hf)\' + \hi-f)\' (0</<oo) (3) 

is the one-sided power spectrum of the detector noise. 
We will show that the bias is most significant near the 
plane of interferometers in the three detector case, agree- 
ing with the observation of reduced resolution in this re- 
gion by |l| . We compute a numerical table of corrections 
for the three-interferometer bias and note that a similar 
procedure may be followed for an arbitrary network of 
detectors. We also find that the bias may be corrected 
through the application of a more detailed parameter es- 
timation [8] , which has the effect of reducing arrival time 
errors. Our application of the least squares estimator 
with bias correction to simulated data is the first test 
of coincident source localization with simulated detector 
noise. 
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Coherent methods also offer promise for source direc- 
tion estimation. The first coherent analysis was described 
by Giirsel and TintoQ. Their approach combines the 
data from a three interferometer network to form a "null 
stream" where the gravitational wave signal should be 
cancelled completely if the assumed source direction is 
correct. The true source direction is estimated by loop- 
ing over a grid of angular locations and minimizing the 
result. Other grid-based approaches to source localiza- 
tion include the maximum likelihood approach of Flana- 
gan and Hughes [31 and the constraint likelihood method 
of Klimcnko et al[J]. The required resolution of the grids 
used in these analyses increases with the maximum fre- 
quency of the search; for LIGO burst searches that extend 
past 1000 Hzfo'] approximately 10^ grid points are neces- 
sary. In practice much coarser grids must be used due to 
computational limitations, potentially leading to missed 
minima and incorrect directional estimates. It should 
also be noted that coherent methods implemented with 
constraints are insensitive to some signal morphologies 
and (small) sky regions. 

The amplitude test described here is similar to the null 
stream approach, but does not require minimization over 
a grid because it works in terms of the root-sum-squared 
strain amplitude 



hrss = y y h{t)^ dt. (4) 

Throughout this paper hrss will refer to the "intrinsic" 
gravitational wave signal at the Earth, i.e. prior to reduc- 
tion by the detector antenna pattern. The test uses the 
hrss values from three non co-located interferometers to 
choose between the two possible source directions given 
by timing considerations and is effective for moderate to 
high SNR. 

The organization of this paper is as follows. Section II 
describes the angular bias in the least squares approach 
to source localization using arrival times from either two 
or three interferometers. In each Monte Carlo sim- 

ulation is pursued and used to characterize and numer- 
ically correct the observed bias. Section III introduces 
the aforementioned amplitude check that can be used to 
differentiate between the two possible source locations 
given by arrival time considerations in the three interfer- 
ometer geometry. Both the time-based source direction 
estimator in (II) and the amplitude check in (III) were 
tested on LIGO- Virgo simulated data and the results arc 
described in Section IV. 



II. Least Squares Source Localization with Bias 
Correction 

In the following we examine the standard least squares 
approach to source localization for both the two and 
three interferometer geometries. We restrict ourselves to 
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FIG. 1: Coordinate system used for source localization. Only 
the polar angle 9 is constrained in the two interferometer 
(IFO) case, while both angles 8 and (j> are constrained in the 
three IFO case. 

studying short duration (< 1 s) "burst" signals so that 
the motion of the Earth may be neglected. We also as- 
sume that the difference in travel time between sites is 
due only to the direction of the source, as explained in 
the introduction. 



A. Two Interferometer Case 

When the network consists of two detectors at different 
locations, the angle 6 between the source unit vector and 
the vector connecting the sites (baseline) is given by 

e = cos-i (^) , (5) 

where r is the time delay, d is the distance between de- 
tectors, and c is the speed of light. This relationship 
is only exact in the absence of noise and constrains the 
source direction to a ring on the sky. When the arrival 
time estimates are affected by noise, the following least 
squares estimator can be used: 

{TT for r < —d/c, 

cos-i(if) for -d/c < T < d/c, (6) 
for T > d/c. 

This is the Maximum Likelihood Estimator and is the 
optimal choice for small arrival time errors. Note that the 
addition of noise introduces timing errors that may lead 
to estimated time delays greater than the light travel time 
between detectors. These unphysical delays are mapped 
to the polar angle that produces the time delay closest 
to what is observed. 

The characteristics of the estimator (6) were studied 
through multiple Monte Carlo simulations. In the first 
simulation we assumed that the distribution of measured 
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Two IFO Bias 




FIG. 2: Least squares angular bias for the Livingston-Hanford 
(Ll-Hl) baseline with a gaussiam arrival time error distribu- 
tion and atH= ""ti =0.1 ms, 1 ms. 



Two IFO Bias witii Multimodal Distribution 




FIG. 3: Least squares angular bias for the Livingston-Hanford 
(Ll-Hl) baseline with a multimodal arrival time error distri- 
bution and crt„= CTtj^ =0.1 ms, 1 ms. 

arrival times for a given source location and detector fol- 
lows a gaussian distribution with mean equal to the true 
arrival time and variance equal to the variance of the 
arrival time estimate. In the second simulation we as- 
sumed that the distribution of measured arrival times 
for a given source location and detector follows a mul- 
timodal distribution consisting of a gaussian main lobe 
and two side lobes each containing half as many points 
as the main lobe. This second scenario is often the case 
for signals with multiple peaks. In both cases we chose 
coordinates as indicated in Figure 1 and used the time 
delay provided by the LIGO Livingston (LI) - LIGO Han- 
ford (HI) baseline. Twenty thousand sets of arrival times 
(with simulated error) from 1800 different sky positions 
were produced via a gaussian random number generator. 
The sky points were spaced by 0.1° on the ^ = 0° merid- 
ian of our symmetric coordinate system. For each sky 



location, the estimator expectation value < 9 > and its 
variance cr? were computed. A nonzero systematic bias 

Bg =<e>-e (7) 

was observed and is plotted against true polar angle 9 
in Figures 2 and 3. The bias in each trial is similar and 
in each case increases with the variance of the arrival 
time distribution. We observed the bias to grow with 
the arrival time uncertainty, and expect it to maintain a 
similar shape for symmetric distributions. The bias is a 
result of two features of the estimator: 

1. ) Mapping of unphysical time delays (|t| > d/c) to 
the two interferometer baseline has the effect of pushing 
the expectation value of the estimated angle towards the 
equator. 

2. ) The nonlinear mapping between the time delay 
r and the angle 9 contributes to the bias through the 
following relation: 

\p{e)d9\ = W)dT\. (8) 

This implies that the distribution p{6) will reflect the 
distribution of time delays p(r) modulated by a sinu- 
soidal Jacobian. Assuming a normal distribution for p(t) 
results in the bias shape seen in Figure 2, while the multi- 
modal distribution discussed above yields the shape seen 
in Figure 3. Note that the bias is a result of the net- 
work geometry and is therefore independent of coordi- 
nate choice. We verified this statement numerically by 
observing the same bias for every sky position when 6 
was chosen along an axis different from the interferome- 
ter baseline. We also verified this property in the three 
interferometer CcLSGj cLS discussed below. 

B. Three Interferometer Case 

In the three interferometer case there are two indepen- 
dent baselines, so a least squares estimator will constrain 
the source direction to two patches on the sky. These 
patches will be mirror images around the plane formed 
by the three detectors and arc indistinguishable because 
they yield identical time delays. The plane of interferom- 
eters is the y-z plane in our chosen coordinate system as 
shown in Figure 1. With these coordinates, we can easily 
determine the output of a least squares estimator given 
three arrival times. The delay times ri, T2 and their un- 
certainties are defined in terms of the three arrival times 
ii, t2, ts: 

Ti=t2-h, (9) 

T2 = t2-ts, (10) 

^-1 = (11) 
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- . /^2 



(12) 



The uncertainties were defined assuming a gaussian 
spread around each true arrival time. Using these defini- 
tions the least squares location estimator (also the max- 
imum likelihood estimator if the arrival times follow a 
gaussian distribution) is equivalent to the following: 

1.) If Inl > di/ck \t2\ > d^/c: 



•£ Tl-dl/c ^ \T2\-d2/c 



r(0,7r) 

(7r,7r) 
(«,^/2) 



j£ T2-d2/c ^ |ri|-di/c 



(13) 



> 



\Ti\~di/c 



2.) If |ril >di/c& \t2\ <d2lc 



(0, tt) if Ti > di/c, 
(tt, tt) if Ti < —di/c 



(14) 



3.) If Inl < di/c & [ral > da/c: 



(^,0) 



(a,7r/2) ifT2>d2/c, 
(tt — a, 37r/2) ifT2<— da/c 



(15) 



4.) If Inl < di/c k (d2/c)(sin(6l)sin(a) + 

COb{0) C0s(q;)) < T2 < c?2/c: 



cos-M^),V2 



(16) 



5.) If Inl < di/c & -d2/c < T2 < 

{d2/c){~ sin{9) sin(a) + cos{6) cos(q:)): 



Least Squares Average Angular Error ( ) 
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FIG. 4: Average angular error in location reconstruction from 
the least squares method applied to the Hl-Ll-Vl network 
when CTti — at2 = (Jts ~ 1 ms. Note that only one hemisphere 
is shown here as we assume the correct sky patch is chosen and 
the result in the other atmosphere will be the mirror image 
of this plot across the y-z plane. 
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cos M ^ ) ,3^2 



(17) 



6.) If Inl < di/c & (d2/c)(sin(6l)sin(a) + 
cos(^) cos(a)) < T2 < (d2/c)(— sin(6') sin(a) -t- 
cos(^) cos(a)): 



cos 



1 / £^ 

di 



cT2/d2 — cos{6) cos(a) 
sin(^) sin(a) 



(18) 



(19) 



The case (1) occurs when both time delays are unphys- 
ical and pins the source location on the baseline that is 
the most standard deviations from physical. The cases 
(2) and (3) occur when one time delay is unphysical and 



FIG. 5: Magnitude of 3 IFO (Hanford (HI), Livingston (LI), 
Virgo (VI)) angular bias when oti — (Jt2 ~ crta ~ 1 ms. This 
magnitude is defined by (22), and the color scale is in degrees. 
Note that only one hemisphere of the sky is shown here, as 
the bias of the other hemisphere is its mirror image across the 
y-z plane. 



places the source location along the appropriate baseline. 
The cases (4) and (5) are a result of analyzing one angle 
first, and place the source at the correct azimuthal an- 
gle if the second time delay is unphysical given the first. 
The last case is the only one where both time delays to- 
gether yield a physical result. Note that the choice of cj) 
is arbitrary when the source is located on a pole. 

The performance of this least squares analysis was eval- 
uated through a Monte Carlo simulation similar to that 
used for the two interferometer case and similar to those 
conducted in [l[. The simulation consisted of 8281 grid 
points spread isotropically across one hemisphere. The 
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FIG. 6: Direction of 3 IFO (HI, LI, VI) angular bias when 
o-ti = crt2 = (Tta = 1 ms. 



hemisphere was chosen to be {tt/2 < 4> < 37r/2, < 6 < 
tt) in our symmetric coordinates (Figure 1); the results 
for the other hemisphere would be a mirror image of 
these. For 3000 iterations at each grid point, the true 
arrival time in each detector was added to normally dis- 
tributed random noise. This gaussian noise was taken to 
have zero mean and variance equal to that of the arrival 
time. 

The average angular error in the location reconstruc- 
tion of this simulation is shown in Figure 4. This error 
was calculated by first computing the angle between the 
true source position (grid point) and estimated location 
in each trial. These angular errors were then averaged 
for each true position (grid point) and turned into the 
map shown in Figure 4. As in the two interferometer 
case, a systematic bias made a significant contribution 
to the overall error. We define the bias in the polar and 
aziniuthal angles as 



uncertainty of 1 ms was assumed for all detectors in both 
plots; in general the bias increases as arrival time uncer- 
tainty increases. 

Note that the location estimation performed by this 
algorithm and the bias correction are independent of co- 
ordinate choice. This is because the estimator places the 
source in the bin that yields time delays closest to what 
is observed regardless of whether the input time delays 
are physical. We verified this property by recomputing 
the bias in an arbitrary coordinate system and observ- 
ing the same results upon rotation back to our original 
coordinates. 

The bias observed in the three interferometer geometry 
can be easily corrected numerically. For a given network 
configuration and set of arrival time uncertainties, one 
can construct a sky map of the effect of the bias. If the 
output of the least squares estimator is then assumed to 
be equal to the expectation value of the estimator, the 
output coordinates can be matched to a term in the bias 
array and the effect of the bias subtracted out. As evi- 
denced by Figures 4 and 5, this correction can reduce the 
mean square errors on the estimated angles significantly. 
The numerical bias correction was applied to LIGO- Virgo 
simulated data and its performance is described in Sec- 
tion IV. 



III. Amplitude Constraints 

Measurements of the gravitational wave signal am- 
plitude at different sites in a network provide valuable 
information for source localization. In the following we 
will describe an amplitude check that resolves the source 
direction degeneracy inherent in the three interferometer 
timing analysis, provided the signal-to-noise ratio is 
sufficiently high in each detector. The antenna patterns 
will be calculated using the method of Thorne I14ll a s 
presented by Yakushin, Klimcnko, and Rakhmanov[10|. 



(20) 



and 



B: =<(!>> 



(21) 



respectively. Note that trials where the source was placed 
on a pole were only counted toward the expectation value 
< ^ >, as 6* is well-defined at these points but is not. 
The total magnitude of the angular bias was calculated 
using all trials and is given by 



M = cos-\P ■ Q), 



(22) 



where P is the unit vector in the direction {< 6 >, < cj) >) 
and Q is the unit vector in the direction of the source. 
A skymap of the value of M for the Hanford-Livingston- 
Virgo (Hl-Ll-Vl) network is plotted in Figure 5, while 
Figure 6 shows the direction of the bias. An arrival time 



A. Detector Response 

In the source frame, the gravitational wave metric h^i, 
can be written in the transverse traceless gauge as 



/O 

h+{t) hy,{t) 

/ix(t) -h+{t) 

Vo 



(23) 



where h^(t) and hx{t) are the two independent polar- 
izations of the source. For gravitational wave bursts of 
short duration, the time dimension can be taken to be 
fixed and we can define the traceless and transverse vec- 
tors 



(24) 
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X 

IFO 

FIG. 7: Coordinate Transformation from IFO to GW frame. 
The IFO frame is defined naturally by the x and y arms of 
the detector. Figure is adopted from [lOl ]. 



and 

/O 1 0\ 
n = 1 0. (25) 

\ oy 

The strain tensor is then 

h{t) = h+{t)ni + hy,{t)n. (26) 

The signal in an interferometer due to an incident gravi- 
tational wave is given by 

S{t)=F+h+{t)+F^h^{t), (27) 

where F+ and Fx are the antenna patterns of the detector 
with respect to the two polarizations. As shown in p^ . 
these coefficients are equal to 

F+ ^ iTr(mR'^mR) (28) 

and 

Fx = iTr(mR'^nR), (29) 

where R is the rotation from the interferometer frame to 
the source frame. This rotation matrix may be decom- 
posed into three Euler rotations, as shown in Figure 7. 
The total rotation is thus 

R = Rz(V'-vr/2)Ry(6l-7r)R^(<?!)), (30) 

where is the polarization angle. Now note that 
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FIG. 8: Efficiency of the two amplitude checks plotted as a 
function of the uncertainty in recovered hrss values. The plus 
check is given by (40) while the minus check is given by (41). 

where h' has the same form as h with 

h'^{t) = cos(2?/')/i+(t) +sin(2V')/ix(0 (32) 

and 

h'^ (t) = - sm(2tl;)h+{t) - cos(2V')/ix (*)■ (33) 

Since the transverse traceless form of h is preserved by 
the final Euler rotation, we may absorb the polarization 
angle into the choice of {h+{t), hx(t)). This choice will 
not affect the response of any interferometer to the grav- 
itational wave source, but does allow us to write new 
antenna patterns as a function of only the source direc- 
tion. 



B. Amplitude Check 

The response of a three interferometer network to grav- 
itational radiation from direction (0, </>) will have the fol- 
lowing form: 

hi{t)=F+sh+{t)+F^^ih^{t) (34) 
h2{t + Ati2) = F+^2h+{t) + Fx,2/ix {t) (35) 

hsit + Atu) = F+^3h+{t) + F^^sh^^it), (36) 

Here the hi terms are the gravitational wave signals, A<i2 
and Ati3 are the time delays from IFO 1 to IFO 2 and 
IFO 3, and F+_i and Fx,i are the plus and cross antenna 
patterns in each detector. Note that the time delays are 
a function of the source direction and we have chosen the 
source coordinate system so as to set the polarization 
angle tjj to tt/2. Combining (34) - (36), we can write one 
signal in terms of the others: 



Rj(^)hR,(^) = h'. 



(31) 



h3{t + Ah3) = A{e,c^)hi{t) + B{e,^)h2{t + Atu), (37) 
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where 



and 



A = 



B ^ 



Fx,lF+^Z — -P+,l-P'x,3 

Fx,iF+^2 — F+^iFx,2 



(38) 



(39) 



Many pipelines produce an estimate of hj-ss (defined by 
(4)) to quantify the strain ampUtude associated with a 
burst event. We can rearrange (37) to include terms of 
this type: 

(40) 

or 

hl,s,3 = iA'+AB)hl,,,, + {B'+AB)hl,^^-ABhl^,^^,_^y 

(41) 

The last terms in (40) and (41) are coherent and refer 
to the sum/ difference of the signals from interferometers 
1 and 2. Both terms are heavily dependent on accurate 
arrival time recovery. The relations (40) and (41) can 
be used to resolve the source direction degeneracy of the 
three detector geometry by comparing the relative er- 
rors of the true location and its mirror image. Either 
test (40) or (41) may be more effective depending on the 
intererometer network and source. A simulation was per- 
formed to determine which check was more efficient for 
the Hl-Ll-Vl network. Linearly polarized gravitational 
wave sources were placed on an isotropic grid of 4060 
points with four different polarization angles chosen ran- 
domly from a uniform distribution between and 7r/2 (Toj 
at each point. Note that the poles and points in the plane 
of the three detectors were excluded as the checks are ir- 
relevant in those directions. The antenna patterns and 
hrss ratios for each trial were calculated and then added 
to noise to see which test held up better as the accuracy 
of the hrss estimates deteriorated. Specifically, each of 
the three hrss values and the sum/difference term were 
multiplied by one plus a normally distributed random 
variable with zero mean and standard deviation equal to 
a specified fraction of the true hrss- This process was 
repeated for each trial and the results were fed to the 
checks (40) and (41) and averaged over all trials to yield 
an efficiency (i.e. the fraction of times the check was suc- 
cessful) for a given relative uncertainty in hrss ■ This was 
repeated for several fractions between and 1, giving the 
plot shown in Figure 8. This plot demonstrates that (40) 
is on average more effective for the chosen network, so 
(40) was chosen for all subsequent tests. 



IV. Application to Simulated Data 

The directional estimator of Section II and the am- 
plitude check of Section III were both applied to LIGO- 
Virgo simulated data. Reference [l3| describes in detail 



how the noise in the three interferometers (HI, LI, VI) 
was produced by filtering stationary Gaussian noise so as 
to resemble the expected spectrum at design sensitivity 
in each detector. Random phase modulation and sinu- 
soids were used to model resonant sources (lines). Six 
types of limearly polarized simulated signals were pro- 
duced using the GravEn algorithm sine gaussians 
with Q=15 and central frequencies of 235 (SG235) and 
820 (SG820) Hz, gaussians of duration 1 (GAUSSl) and 
4 (GAUSS4) ms, and two families of Dimmelmeier-Font- 
Mueller (DFM) supernova core collapse waveforms [Tl| 
with parameters A=l, B=2, G=l and A=2, B=4, G=l. 
These waveforms are shown in Figure 9. The injections 
were configured as above to produce an isotropic distri- 
bution of 3960 points on the sky, with each of these grid 
points being tested with 4 different polarization angles 
(-0 = 0, 7r/8, 7r/4, 37r/8). The waveforms were delayed 
and scaled appropriately and added to the noise of each 
detector for analysis. The process was repeated for each 
of the six waveforms. 

A stand-alone parameter estimation Q was used to 
determine the arrival times, their uncertainties, and the 
hrss values in each trial. The parameter estimation pro- 
cessed input data by first applying high pass (100 Hz) 
and linear predictor (whitening) zero-phase filters. The 
noise power spectrum was determined by averaging the 
spectra of an interval before and an interval after the in- 
jected signal. This power spectrum was subtracted from 
the spectrum of the interval containing the signal and the 
result divided by the whitening filter frequency response. 
The measured amplitude /i™^'^"'' was computed over the 
band from 100 to 1000 Hz as 



•Anhand 



1000//2 



(42) 



subject to a threshold on the signal-to-noise ratio. Specif- 
ically, the data was broken into N frequency bins and 
the content in a given bin [/i-i, fi] {i = 1..N) was only 
counted if 




WW 
s{f) 



df>9. 



(43) 



This SNR threshold was chosen to emphasize compo- 
nents with moderate visibility above the noise floor while 
reducing spurious excess power. The estimated signal 
power spectrum was also used to determine the frequency 
band that contained the middle 90% of the signal power. 
The original whitened time series was zero-phase band- 
passed according to this estimate and its maximum taken 
as the arrival time. The arrival time uncertainty was 
taken to be the Cramer-Rao Lower Bound of the arrival 
time estimator. 

The directional estimator consisted of the standard 
least squares estimator with bias correction as described 
in Section II. The bias correction array had directional 
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FIG. 9: Simulated waveforms used for amplitude check test- 
ing. 
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io-i« 


SG235 


0.5415 


0.8924 


0.9470 


0.9947 


0.9994 


SG820 


0.0252 


0.6438 


0.8076 


0.9803 


0.9981 


GAUSSl 


0.3210 


0.8215 


0.9071 


0.9904 


0.9993 


GAUSS4 


0.5537 


0.8936 


0.9458 


0.9949 


0.9995 


DFM A1B2G1 


0.0160 


0.6416 


0.8097 


0.9799 


0.9979 


DFM A2B4G1 


0.3338 


0.8238 


0.9097 


0.9910 


0.9992 



TABLE I: Fraction of events with sufficient energy for pa- 
rameter estimation to estimate parameters. Note tfiat this 
fraction is larger than what would be detected by a search 
algorithm tuned to the same frequency range [l3] due to in- 
clusion of low energy events that would fall below standard 
detection thresholds. The performance of the directional esti- 
mator presented here is therefore a conservative estimate due 
to the inclusion of these low energy events. 

resolution of 2 degrees azimuthally and 90 bins evenly 
spaced in the cosine of the polar angle. Its timing reso- 
lution was 0.1 ms in the arrival time uncertainty at each 
detector. Trials were conducted with intrinsic inband 
hrss varying from lO^^^Hz^^/'^ to 10~^^Hz~^^^, corre- 
sponding to realistic to loud astrophysical sources. These 
amplitudes were scaled down by the antenna pattern be- 
fore being added to the detector noise, on average mul- 
tiplying the incident hrss by 0.38. The plot in Figure 10 
shows the average angular error for all waveforms as a 
function of intrinsic inband hrss ■ This error is defined as 

A = cos-\V ■ R), (44) 

where V is the unit vector in the direction of the es- 
timate and R is the unit vector in the direction of the 
source. The average was taken over all sky locations and 
polarization angles where a signal was visible to the pa- 
rameter estimation in all three detectors, and was found 
to be virtually identical with and without the bias cor- 
rection. As can be seen from Figure 11, a few outlying 
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FIG. 10: Performance of least squares angular estimator on 
LIGO- Virgo Project lb simulated waveforms for varying in- 
trinsic inband hrsa- In all cases the effect of the bias correction 
was negligible due to sufficiently small arrival time uncertain- 
ties. 

samples caused the average error in each case to be signif- 
icantly higher than most individual errors. These outliers 
were due to incorrect peak time estimates from the pa- 
rameter estimation. In most cases the algorithm picked a 
secondary peak of the signal (see Figure 12), causing the 
arrival time estimate to be off by a few milliseconds and 
resulting in a large error in the directional estimate. The 
SG820 and DFM A1B2G1 waveforms were particularly 
susceptible to this, as shot noise becomes firmly estab- 
lished at the higher frequencies where these signals are 
concentrated. The lack of impact of the bias correction 
can be attributed to the relatively small timing uncer- 
tainties (approximately 0.1 ms) given by the parameter 
estimation. The bias is negligible for such small timing 
errors, so the correction did not lead to substantial im- 
provement in the trials undertaken. Even with weaker 
input signals, the bias correction was seen to be inconse- 
quential as the arrival time uncertainties remained low. 
It should also be noted that the errors in the arrival time 
estimates were observed to follow multi-modal distribu- 
tions in waveforms with multiple peaks, as the most sig- 
nificant errors were a result of the parameter estimation 
choosing the wrong peak. 

The observed angular uncertainties can be placed in 
context through comparison with the classical limit com- 
monly used in radio astronomy [l^ [l^, [l3|. The classical 
limit 9 is given by 



where A is the wavelength of the radiation being con- 
sidered and D is the baseline of the network. Using the 
longest such baseline in our network (Hl-Vl) and assum- 
ing a frequency in the middle of our sensitive band (550 
Hz) gives a classical limit for our study of Oqw ^ 3.8°. 
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FIG. 11: Histogram of angular errors for all waveforms with 
an intrinsic inband hraa of 10~^'^ Hz~^^^ . Note the logarithmic 
scale of the y-axis. 
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FIG. 12: Example of difficulty in arrival time estimation. The 
parameter estimation picks the largest peak in the whitened, 
bandpassed data but this may not be the primary peak due 
to noise. 

This differs by less than an order of magnitude from the 
observed average angular error for all waveforms at an 
intrinsic inband hrss of 10~^°Hz^^^, which corresponds 
to a loud astrophysical source (SNR ~ 20). 

The same distribution of points, polarization angles, 
and signal strengths was used to determine the efficiency 
of the amplitude check (40) as a function of input sig- 
nal amplitude. The simulations were again generated in 
terms of intrinsic inband hrss and scaled down by their 
antenna factors before being added to the data. Note 
that the reported efficiencies are an average over all di- 
rections and polarization angles where the signal was 
detected in all three interferometers. Two sets of tri- 
als were conducted: one where the exact arrival times 
were assumed to be known for forming the last (coher- 
ent) term in (40) and one where the parameter estimation 
peak times were used for its determination. The efficien- 
cies were only slightly better with the former "aligned" 
approach, showing that the peak times from the pa- 
rameter estimation are sufficient to accurately calculate 
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FIG. 13: Average efficiency of amplitude check (40) as a func- 
tion of intrinsic in-band ftrss for each of the sample waveforms. 
The hras values shown here were scaled down further by the 
antenna factors before entering the data. Note that these 
curves were generated assuming that the time delays between 
IFOs 1 and 2 were known exactly. 
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FIG. 14: Average efficiency of amplitude check (40) as above 
except with the last term of (40) being determined using pa- 
rameter estimation peak times. 

the coherent term. The efficiencies for the SG820 and 
DFM A1B2G1 waveforms were slighly below those of the 
other waveforms due to their power being concentrated 
at higher frequencies where the interferometer noise is 
worse. 



V. Conclusions 

This work focused on two aspects of coincident source 

localization. The first was a least squares analysis that 
used signal arrival times and their uncertainties to es- 
timate the direction of the source. This estimator was 
studied and found to be biased in both the two and three 
interferometer cases, with the effect increasing with ar- 
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rival time uncertainty. The cause of the bias was seen 
to be twofold as both the allocation of unphysical time 
delays to their most likely direction and the nonlinear 
mapping between the time delays and source coordinates 
contributed to the effect. The bias will exist for networks 
with an arbitrary number of detectors, and can be de- 
termined numerically through Monte Carlo simulations 
with a minimization routine such as that described 
in [l|. In each case the bias can be corrected numeri- 
cally with the correction being relevant for arrival time 
uncertainties greater than about 1 ms. Studies on real 
data "T] indicate that such uncertainties are probable 
using conventional detection algorithms. Our analysis in- 
dicates that the bias may also be corrected through the 
reduction of arrival time errors provided by a thorough 
parameter estimation such as the one used hereQ. 

The second piece of coincident analysis described here 
was an amplitude check, applicable to both linear and 
multiply polarized signals, that can be used to resolve the 
source location degeneracy inherent in the three interfer- 
ometer geometry. While this analysis is not coincident in 
a strict sense (one term requires combined data streams), 
it is extremely simple and so provides a lightweight alter- 
native to more sophisticated methods. As shown in Sec- 
tion IV, the check is effective for realistic to loud sources, 
particularly those where the peak time may be recovered 
accurately. 

One potential application for these source localization 
techniques would be the generation of skymaps of back- 
ground and foreground triggers. In this case the ampli- 
tude check would improve the efficiency of a distribu- 



tional test based on sky locations even for cases where 
its efficiency is only slightly better than 50%. We would 
expect such maps to be isotropic based on current detec- 
tor sensitivity (i.e. non-gravitational wave triggers dom- 
inate) but in the the future such plots may give us a 
picture of the gravitational wave sky. 

Methods such as those described above are useful in 
that they are inexpensive computationally and quite ac- 
curate in most cases. However they may struggle if a 
source has similar amplitudes in both the -t- and x po- 
larizations as well as a delay between the peaks of the 
waveforms. Due to the differing antenna responses of 
the members of an interferometer network, peaks corre- 
sponding to the different polarizations may be recorded 
in different detectors, leading to an inconsistency in ar- 
rival time estimates that cannot be accounted for with 
coincident methods. The likelihood of this situation is 
currently under investigation, as is the performance of 
the algorithms presented here on real data with randomly 
polarized signals. 



VI. Acknowledgements 

LIGO was constructed by the California Institute 
of Technology and Massachusetts Institute of Technol- 
ogy with funding from the National Science Founda- 
tion and operates under cooperative agreement PHY- 
0107417. This paper has LIGO Document Number 
LIGO-P080003-00-Z. 



[1] F. Cavalier et al, Phys. Rev. D 74, 082004 (2006). 
[2] Y. Giirsel and M. Tinto; Phys. Rev. D 40, 3884 (1989). 
[3] E. E. Flanagan, S. A. Hughes, Phys. Rev. D 57, 4535 
(1998) 

[4] S. Klimenko, S. Mohanty, M. Rakhmanov, G. Mitsel- 
makher, Phys. Rev. D 72, 122002 (2005)^ 

[5] Chatterji et al. Class. Quantum Grav. gr-qc/0605002| 

[6] I. Yakushin: Search for Gravitational Wave Bursts in 
LIGO's S5 Run, LIGO-G060170 (2006)*. 

[7] N. Arnaud et al, Phys. Rev. D 67, 062004 (2003). 

[8] M. Zanolin et al., LIGO-T060226-00-R (2008)*. 

[9] Abbot B et al (LSC) 2006 Class Quantum Grav. 23 S29- 
39 

[10] I. Yakushin, S. Klimenko, and M. Rakhmanov: MDC 
Frames for S2 Burst Analysis, LIGO-T040042-00-Z-00- 
D (2004)*. 

[11] H . Dimmelmeier, J. Font and E. Mueller, Astron. Astro- 

phys. 393 , 523-542 (2002). 
[12] F. Beauville et al. 2008 Class. Quantum Grav. 25 045002 

(32pp). 

[13] A. Stuver: GravEn Simulation Primer, LIGO-T040020- 
01-Z (2004)*. 



[14] 

[15] 
[16] 
[17] 

[18] 
[19] 



K.S. Thorne, in 300 Years of Gravitation, edited by S. 
Hawking and W. Israel (Cambridge Universtiy Press, 
Cambridge, 1987), chap. 9, pp.330-458. 
Carroll and Ostlie: An Introduction to Modem Astro- 
physics (Addison- Wesley, Reading MA, 1996). 
J. Kraus: Radio Astronomy, Second Edition (Cygnus- 
Quasar Books, Powell OH, 1986). 

Thompson, Moran, and Swenson: Interferometry and 
Synthesis in Radio Astronomy, Second Edition (John Wi- 
ley & Sons, Inc., New York, 2001). 

F. Beauville et al. 2008 Class. Quantum Grav. 25 045001 
(22pp). 

The polarization angles were chosen in this manner to 
avoid redundancy. The gravitational wave oscillations 
from -i/) = to ^/) = TT are identical to those from tp = n 
to i/> = 27r, while only a minus sign differentiates the an- 
tenna patterns from ijj = Q to ijj = -k /2 and t/j = 7r/2 to 
ijj = TT. The minus sign should not affect the hrss values 
used in the check, so only polarizations between i/) = 
and ip = n/2 were considered. 



* Available at http://docuserv.ligo.caltech.edu/ 



